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D.  P-  O’Leary 
G.  W.  Stewart 


1.  Introduction 

This  is  a  summary  of  work  accomplished  under  Grant  AFOSR-82-0078.  The 
purpose  of  this  effort  is  to  develop  realistic  algorithms  for  matrix  computations 
on  parallel  computers.  It  has  been  long  observed  that  the  usual  algorithms  of 
numerical  linear  algebra  contain  a  great  deal  of  inherent  parallelism.  For  exam¬ 
ple,  if  the  arithmetic  operations  that  can  be  performed  in  parallel  in  Gaussian 
elimination  are  actually  so  executed,  the  time  to  decompose  an  nXn  matrix  is 
reduced  from  order  t?  to  n  Only  recently,  with  the  emergence  of  cheap,  small 
microcomputers,  has  it  become  feasible  to  exploit  this  parallelism  on  anything 
but  a  trivial  scale. 

At  the  Department  of  Computer  Science  at  the  University  of  Maryland, 
there  is  under  development  a  parallel  system,  called  the  ZMOB,  consisting  of 
256  micro-processors  connected  on  a  conveyor  belt.  This  belt  is  so  fast  and  its 
architecture  b  such  that  any  two  processors  can  communicate  without  interfering 
with  the  communications  of  other  pairs  of  processors.  Thus  the  ZMOB  is  an 
ideal  tool  for  simulating  an  arbitrarily  connected  network  of  computers. 

This  feature  of  the  ZMOB  is  particularly  useful  in  investigating  parallel 
matrix  algorithms.  As  was  noted  above,  there  is  much  parallelism  in  most  current 
matrix  algorithms.  However,  to  exploit  it,  information  must  be  moved  from  pro¬ 
cessor  to  processor.  This  constitutes  the  chief  bottleneck  in  parallel  matrix  algo¬ 
rithms;  interconnections  between  processors  are  expensive,  and  in  a  practical  sys¬ 
tem  one  can  assume  only  a  limited  amount  of  connectivity.  The  ZMOB  provides 
a  means  of  testing  and  comparing  different  types  of  interconnections,  since  all  one 
has  to  do  b  not  use  the  rich  connections  provided  by  the  ZMOB  conveyor  belt. 
Thus  our  proposal  b  to  use  the  ZMOB  to  design  and  test  networks  for  parallel 
matrix  computations. 


2.  Progreu  to  Date 

Our  research  b  proceeding  in  three  stages.  First,  decide  on  a  suitable  way  of 
connecting  and  synchronizing  processors  for  parallel  matrix  computations. 
Second,  design  and  build  a  communications  system  to  realize  thb  network  on  the 
ZMOB.  Third,  code  matrix  algorithms  for  the  system,  and  experiment  with 
them.  In  addition,  we  must  install  and  test  the  floating-point  processors  which 
were  requested  as  part  of  the  initial  grant  period.  In  this  section  we  shall  take  up 
each  of  these  points  in  turn.  _ _  ,  ,  _ -Tvr-  '  ^ 
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A  greater  part  of  our  research  coucerns  two  dimensional  arrays  of  processors, 
and  we  have  made  considerable  progress  in  this  area.  It  is  highly  desirable  to  be 
able  to  restrict  the  connections  in  such  an  array  to  lines  between  adjacent  proces¬ 
sors,  since  this  is  the  simplest  and  most  easily  implemented  of  networks.  We 
have  observed  that  not  only  can  many  matrix  algorithms  be  implemented  on  such 
a  network,  but  also  the  processors  can  be  synchronized  by  the  flow  of  data  in  the 
network,  without  any  need  of  outside  control.  A  sketch  of  how  such  networks 
operate  was  given  in  the  first  renewal  proposal.  Here  we  just  list  some  of  the 
advantages  of  the  approach. 


1.  The  interconnections  are  simple  and  realizable. 

2.  Each  processor  can  operate  asynchronously. 

3.  The  same  program  can  be  used  on  each  processor. 

4.  Many  matrix  algorithms  fit  naturally  into  this  scheme. 

5.  The  approach  provides  a  natural  way  of  dealing  with  array  overflow; 
i.e.,  the  case  where  the  size  of  the  matrix  exceeds  the  size  of  the  array 
processors. 

In  order  to  support  the  network,  we  are  building  a  comunications  system  to 
pass  information  from  processor  to  processor.  This  system  will  be  invoked  from  a 
high  level  programming  language,  and  it  will  permit  multi-processing  on  a  single 
processor.  This  latter  feature  is  necessary  to  cope  with  array  overflow.  The  core 
of  the  operating  system  has  been  programmed  and  has  been  used  to  perform 
small  matrix  computations  on  the  ZMOB.  We  shall  begin  testing  code  previously 
written  for  the  system. 

Although  much  of  our  current  effort  is  devoted  to  building  a  system  for  test¬ 
ing  parallel  matrix  algorithms,  we  are  abo  designing  new  parallel  algorithms  for 
important  matrix  processes.  In  particular  we  have  developed  a  promising  algo¬ 
rithm  for  the  solution  of  the  non-Hermitian  eigenvalue  problem.  The  method  is 
based  on  a  Jacobi-like  iteration  to  reduce  a  matrix  to  upper  triangular  form  by 
unitary  transformations.  It  is  numerically  stable  and  parallelizes  readily.  Prelim¬ 
inary  experiments  indicate  that  it  will  be  effective  for  a  wide  class  of  eigenvalue 
problems. 

Work  has  also  been  done  on  parallel  algorithms  for  solving  sparse  matrix 
problems  that  have  a  sparsity  structure  corresponding  to  a  grid  of  points  con¬ 
nected  to  up  to  eight  nearest  neighbors.  Such  problems  arise  in  discretization  of 
elliptic  partial  differential  equations,  network  problems,  and  image  processing. 
Various  three-colorings  of  the  graph  and  corresponding  numberings  of  mesh 
points  have  been  devised  so  that  an  iteration  of  a  relaxation  algorithm  such  as 
Gauss-Seidcl  or  SOR  can  be  executed  with  Parallelism  comparable  to  the  Jacobi 
algorithm,  without  degradation  of  convergence  rate. 
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We  have  also  been  investigating  algorithms  for  determining  the  equilibrium 
vector  of  nearly  uncoupled  Markov  chains.  These  chains  arise  naturally  in  the 
stochastic  modeling  of  computer  systems.  We  have  analysed  the  properties  of  a 
highly  parallelizable  method  based  on  a  combination  of  aggregation  and  the  block 
Gauss'Seidel  method. 

Finally,  we  have  coded  a  test  package  for  the  floating  point  processors,  so 
that  they  may  be  quickly  incorporated  into  the  individual  ZMOB  boards. 

Although  at  this  time  the  project  is  still  in  a  developmental  state,  we  have 
given  several  talks  on  our  work  and  have  prepared  papers  for  publication.  These 
are  listed  in  Appendix  A. 
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Appendix  A 


I.  Technical  Reports 


(1)  G.  W.  Stewart,  Computing  the  CS  Decompoeition  of  a  Partitioned  Orthonor' 
mat  Matrix,  TR-1159,  May,  1082. 

This  paper  describes  an  algorithm  for  simultaneously  diagonalizing  by 
orthogonal  transformation  the  blocks  of  a  partitioned  matrix  having  ortho¬ 
normal  columns. 

(2)  G.  W.  Stewart  A  Note  on  Complex  Division,  TR-1206,  August,  1982. 

An  algorithm  (Smith,  1962)  for  computing  the  quotient  of  two  complex 
numbers  is  modified  to  make  it  more  robust  in  the  presence  of  underflows. 

(3)  D.  P.  O’Leary,  Solving  Sparse  Matrix  Problems  on  Parallel  Computers,  TR- 
1234,  December,  1982. 

This  paper  has  a  dual  character.  The  flrst  part  is  a  survey  of  some  issues 
and  ideas  for  sparse  matrix  computation  on  parallel  processing  machines.  In 
the  second  part,  some  new  results  are  presented  concerning  efficient  parallel 
iterative  algorithms  for  solving  mesh  problems  which  arise  in  network  prob¬ 
lems,  image  processing,  and  discretization  of  partial  differential  equations. 

(4)  G.  W.  Stewart,  A  Jaeobi-like  Algorithm  for  Computing  the  Sehur  Decomposi¬ 
tion  of  a  Non-IIermitian  Matrix,  TR-1321,  August,  1083. 

This  paper  describes  an  iterative  method  for  reducing  a  general  matrix  to 
upper  triangular  form  by  unitary  similarity  transformations.  The  method  is 
similar  to  Jacobi's  method  for  the  symmetric  eigenvalue  problem  in  that  it 
uses  plane  rotation;  to  annihilate  off-diagonal  elements,  and  when  the  matrix 
b  Hermitian  it  reduces  to  a  variant  of  Jacobi’s  method.  Although  the 
method  cannot  compete  with  the  QR  algorithm  in  serial  implementation,  it 
admits  of  a  parallel  implementation  in  which  a  double  sweep  of  the  matrix 
can  be  done  in  time  proportional  to  the  order  of  the  matrix. 


II.  Technical  reports  in  preparation 
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(1)  D.  P.  O’Leary  and  G.  W.  Stewart,  Data  Flow  Algorithma  for  Matrix  Compv- 
tationa, 

(2)  G.  W.  Stewart  and  R.  van  de  Geijn,  VMOB:  Virtual  ZMOB, 

(3)  D.  McAllister,  G.  W.  Stewart,  and  W.  J.  Stewart,  A  Two-Stage  Algorithm 
for  Nearly  Uncoupled  Markov  Chains, 

(4)  D.  P.  O’Leary,  Block  Preconditioninga  for  Parallel  Computations, 

III.  Presentations  during  1083 

(1)  D.  P.  O’Leary,  Sohing  Mesh  Problems  on  Parallel  Computers, 

Bel!  Laboratory,  Murray  Hill,  N.J.,  January,  1983 

IBM  T.  J.  VV’atson  Laboratory,  Yorktown  Heights,  N.Y.,  January,  1983. 

(2)  G.  W.  Sewart,  A  Jacobi-like  Algorithm  for  Computing  the  Schur  Decomposi¬ 
tion  of  a  Non-Hermitian  Matrix  (invited).  Symposium  on  Numerical  Analysis 
and  Computational  Complex  Analysis,  Zurich,  Switzerland,  August,  1983. 

(3)  G.  W.  Stewart,  The  Structure  of  Nearly  Uncoupled  Markov  Chains  (invited). 
International  Workshop  on  Systems  Modeling,  Pisa,  Italy,  September,  1083. 

(4)  G.  W.  Stewart,  Data  Flow  Algorithms  for  Paradlel  Matrix  Computation's 
(invited),  SL4M  Conference  on  Parallel  Processing  for  Scientific  Computing, 
Norfolk,  VA,  November,  1983. 

(5)  D.  P.  O’Leary,  Parallel  Compulations  for  Sparse  Linear  Systems  (minii^^ympo- 
sium  invitation),  SIAM  1983  Fall  Meeting,  Norfolk,  VA,  November,  19^. 

(6)  D.  C.  Fbher,  Numerical  Computations  on  Multiprocessors  with  Only  Local 
Communications  (poster  session),  SIAM  Conference  on  Parallel  Processing 
for  Scientific  Computing,  Norfolk,  VA,  November,  1983. 
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